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Aerodynamic shape optimization was performed on an isolated axisymmetric plug nozzle 
sized for a supersonic business jet. The dual-stream concept was tailored to attenuate 
nearfield pressure disturbances without compromising nozzle performance. Adjoint-based 
anisotropic mesh refinement was applied to resolve nearfield compression and expansion 
features in the baseline viscous grid. Deformed versions of the adapted grid were used for 
subsequent adjoint-driven shape optimization. For design, a nonlinear gradient-based 
optimizer was coupled to the discrete adjoint formulation of the Reynolds -averaged Navier- 
Stokes equations. All nozzle surfaces were parameterized using 3 rd order B-spline 
interpolants and perturbed axisymmetrically via free-form deformation. Geometry 
deformations were performed using 20 design variables shared between the outer cowl, 
shroud and centerbody nozzle surfaces. Interior volume grid deformation during design was 
accomplished using linear elastic mesh morphing. The nozzle optimization was performed at 
a design cruise speed of Mach 1.6, assuming core and bypass pressure ratios of 6.19 and 3.24, 
respectively. Ambient flight conditions at design were commensurate with 45,000-ft standard 
day atmosphere. 


I. Introduction 

O verland noise is a challenge high speed aircraft must overcome to gain regulatory permission to operate in the 
supersonic regime. Recent conceptual design studies have shown moderate reductions in perceived noise level 
can be achieved through careful aerodynamic tailoring of the airframe. 1-4 These studies often neglect powered 
propulsion, discounting the presence of exhaust jets in nearfield pressure signatures and propagated sonic boom 
metrics. Furthermore, on-going scale experiments 5,6 have validated the sonic boom performance of promising low- 
boom aircraft candidates, but also lack powered nacelles and provide no insight into propulsion- airframe interaction 
effects. While the amplitudes of discrete pressure disturbances from the propulsion system have been shown minor 
relative to those induced by the airframe 7 , experimental and computational studies 8 " 10 both suggest shock and plume 
interactions become less negligible for vehicles designed to meet a low-boom standard. In particular, introducing 
propulsion effects into an optimized airframe boom signature has been found in multiple cases to compromise the 
low-boom requirement. 7,10 

To minimize powered propulsion contributions to the overall aircraft boom signature, aerodynamic shape 
optimization has been applied to an axisymmetric dual- stream supersonic nozzle. The objective of this research is to 
mitigate nearfield pressure waveforms produced by the isolated nozzle jet without compromising nozzle 
performance. This successful capability will simplify propulsion- airframe integration for ultra low-boom aircraft and 
enable propulsion system boom contributions to be further minimized in the presence of an airframe. 

While this analysis does not directly consider perceived noise level by a ground observer, the nearfield pressure 
signature is a known primary driver for sonic- wave propagation theory. Signal attenuation attributed to atmospheric 
nonlinearities, thermo-viscous absorption and molecular relaxation are acknowledged but considered beyond the 
scope of this study. 

Aerodynamic analysis and design optimization were performed using the computational fluid dynamics (CFD) 
package FUN3D 11 . FUN3D is a three-dimensional, node-centered, Euler and Navier-Stokes flow solver for fully 
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unstructured grids. Tailored for flow analysis and design, the code offers a discretely consistent adjoint methodology 
enabling sensitivity calculations necessary for gradient-based mesh adaptation and design optimization. 

II. Problem Definition 

A. Baseline Geometry 

Figure 1 is an isometric view of the baseline geometry, comprised of a centerbody, inner shroud and outer cowl. 
The geometry was adapted from an axisymmetric expansion plug concept 10 developed by the Gulfstream Aerospace 
Corporation (GAC) with dual-annular high-flow nacelle bypass exhaust 12 ’ 13 . The original nozzle was created using 
1-D gas dynamics and method-of-characteristics (MOC) combined with streamtube shaping, an inverse-design 
approach for tailoring extensions to the nozzle expansion surfaces to attenuate pressure field disturbances. 10 The 
core stream passes between the centerbody and inner shroud, while the pressurized high-flow bypass stream, which 
serves as an aero-acoustic shield, flows between the shroud and outer cowl. During operation, a mechanical actuator 
linearly translates the centerbody to adjust the effective throat area and maintain performance at off-design. This 
study only considers nozzle shape optimization of the cruise flight configuration. 

The original nozzle reported by GAC contained a centerbody extended significantly longer than that shown in 
Fig. 1. An extended plug offers the benefits of a reduced momentum wake in the plume and less flow turning near 
the plug tip. Despite these aerodynamic gains, the added weight penalty and viscous losses along the centerbody 
surface were reasons a shortened plug variant was created for this study. 



Fig. 1. Isometric view of the baseline plug nozzle. 


B. Boundary Conditions 

Despite axisymmetric geometry, current limitations of the FUN3D design optimization framework required 
analysis of a three-dimensional quarter- symmetry sector. The 90-degree computational domain in Fig. 2 was 
constructed with y- and z-planar symmetry prescribed on the vertical and horizontal surfaces of the domain. The x- 
axis corresponds to the axial flow direction and is aligned with the nozzle centerline. A farfield Riemann invariant 
condition was assigned to the outer cylindrical boundary. Subsonic inflow was specified with total temperature and 
total pressure ratios relative to freestream static conditions for the core and bypass streams. A freestream condition 
was prescribed on the primary inflow face and a static backpressure permitting supersonic flow that extrapolated 
from the interior of the domain was applied at the downstream outflow boundary. All nozzle surfaces were modeled 
as viscous walls with strong enforcement of no- slip. 

Bypass and core stream boundary conditions were supplied based on turbine exit conditions provided by GAC. 
The core nozzle pressure ratio (NPR) or ratio of the total pressure to freestream static pressure was 6.19. The core 
temperature ratio (NTR) or total core temperature normalized by freestream static temperature was 2.43. The bypass 
stream NPR was prescribed as 3.24 with corresponding NTR set at 1.51. The freestream reference conditions for 
cruise corresponded with 45,000-ft standard day atmosphere and a cruise speed of Mach 1.6. The Reynolds number 
per unit grid length used for solver non-dimensionalization was 7.8 x 10 6 . Flowfield conditions are summarized in 
Table 1 where M^, P^ and T^ represent freestream Mach number, static pressure and temperature, respectively. For 
reference, the baseline nozzle produces approximately 88,000-N (19,800-lbf) thrust at cruise. 
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TABLE 1 : Flowfield Conditions 


Parameter 

Value 

Units 

M ro 

1.6 

- 

Tco 

216.65 

K 

Poo 

14747.89 

Pa 

R^unit 

7.8 x 10 6 

- 

NPRcore 

6.19 

- 

NTR core 

2.43 

- 

NPRbypasS 

3.24 

- 

NTRbypass 

1.51 

- 



Farfield 



Supersonic 

Outflow 


C. Optimization Problem Formulation 

The formal optimization problem may be stated as: 


Minimize: 



( 1 ) 
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Where / is the objective function defined as the sum of squared differences between the local static pressure 

i i * 

ratio — obtained along a linear signal at point i and a pre-defined target static pressure ratio — . For low-boom 

Poo Ij Poo '[ 

aircraft shape optimization, a nearfield pressure target is typically specified from a desired low amplitude ground 
signature based on wave propagation mechanics. The ground signature is defined from an equivalent area target 
distribution compatible with the vehicle configuration. The characteristic length of the target signature pre-defines 
the size of the body generating the nearfield pressure disturbance. Unlike vehicle design where aircraft length is pre- 
determined by mission requirements, the geometric nozzle length is unknown a priori. Rather than identify a 
nearfield pressure target to match using the equivalent area approach, a zero disturbance local pressure field was 
defined over the entire target distribution. While this target distribution is non-physical, it places no assumed 
constraints on nozzle length. 

The nozzle geometry was constrained in parametric space to maintain core and bypass stream throat areas and 
avoid compromising nozzle performance during optimization. Thrust was computed using the propulsion 
performance package in FUN3D that was recently upgraded and verified against the WIND 14,15 structured flow 
solver. Thrust was determined by application of the momentum theorem to a control volume about the nozzle. 
Viscous skin friction and pressure forces were integrated over each solid nozzle surface boundary while axial 
pressure and momentum fluxes were computed for the core and bypass stream inflow boundaries. Summation of the 
force and flux contributions combined to give total thrust for the isolated nozzle. If integrated with an airframe, only 
throttle dependent surfaces would by bookkept with engine thrust. For example, the viscous and pressure forces on 
the outer cowl surface would be excluded from the engine thrust computation and included as an aircraft drag 
penalty. Given this study only considers an isolated nozzle, the forces acting on the outer cowl must be included 
when computing overall nozzle performance. 

Axial pressure forces and fluxes within FUN3D were internally normalized by the freestream reference Mach 
number, pressure and speed of sound as shown in Eq. 2, where F c is the normalized thrust component for each 
surface or inflow boundary. The ratio of specific heats is denoted as y. 

Thrust = F c (2) 

III. Geometry Parameterization 

To support this work, a new parameterization was developed specifically for axisymmetric geometries and 
integrated with FUN3D. The implementation used free-form deformation (FFD) with 3 rd order B-splines. FFD was 
selected to parameterize changes to the initial geometry, rather than the geometry directly, resulting in a reduced 
number of degrees of freedom and simplified design space. The selection of B-splines provided a straightforward 
means for computing analytic sensitivities. 

The FFD approach begins with a B-spline interpolant, defined as a function including parametric coordinate, u , 
and control point vector, C, of length n , given by Eq. 3. In the definition, b Jp (u) is a polynomial basis function of 
degree p with knot vector t , and defined recursively by Eq. 4. It should be noted that knot vector selection is non- 
arbitrary and controls overall spline behavior. 16 Knot spacing was distributed uniformly for each control spline. 

n 

spline (u, C) = ^ Cjb jp (u ) (3) 

7=0 

if tj <u < tj +1 and tj < t j+1 
otherwise 

t j L , t j+p+1 -U ^ ^ (4) 

_ f "f . _ . bj+i t p-i&) 

l j l j+p + i l j+i 

The axi-symmetric assumption permits application of 1 -dimensional B-splines, which means that u is a scalar, 
while C is a vector of 2-dimensional coordinates. C is defined here using polar coordinates with r (radial) and x 
(axial) components. For each individually parameterized nozzle body (centerbody, cowl and shroud), axial locations 
of the initial and final control points were required to correspond with the axial locations of initial and final points of 
the discrete baseline geometry. 

Once the B-spline interpolant was defined, a mapping was found to link the parametric space of the B-spline 
with the physical space of the discrete geometry. Each discrete point on the geometry, P u was mapped to the B- 
spline parametric coordinate, u h by solving the non-linear system in Eq. 5. 

4 

American Institute of Aeronautics and Astronautics 


bj iP (u) 


h; 


hv = 


(5) 


spline (u{) x = P ix 


dP- 

Given the mapping from P t to u h the total geometric derivatives, were computed using Eqs. 6 and 7. Notably, 


Eq. 5 is constant with respect to C. 


dspline 

dCj 


n 

^ bj, v (u) 

i = 0 


( 6 ) 


dP { 

dCj 


dP t dspline 


dspline dCj 


Ui 


n 

y, bj >p (Ui) 

i = 0 


(7) 


B-spline interpolants were applied in two different ways to individually parameterize the nozzle surfaces. The 
nozzle centerbody, a closed solid, was deformed using a single B-spline to control radial and axial deformations 
using Eqs. 8 and 9. 


R t = R i>0 + R re f \spline(u it Cj)] (8) 

X t = spline(u if X ref Cj ) (9) 

R i>0 denotes the radial coordinate of the non-deformed geometry at discrete point i. R re f and X ref are parameters used 
to scale the magnitude of deformations in physical space relative to deformations of the spline. The size of the 
resulting design space is a function of the number of control points in C. For the centerbody parameterization, each 
control point contributes two variables (one radial and one axial degrees of freedom) to the design space. 

Application of Eq. 7 yields the following derivatives of the deformations in the axial and radial directions with 
respect to C: 


dR t dspline 

dq = Rref dCj 


dX t dspline 

dq =Xref dCj 


( 10 ) 


(ii) 


The cowl and shroud surfaces, hollow shell bodies, were each divided into inner and outer surfaces and 
parameterized using a pair of independent B- splines. The first spline was used to control shell centerline profile and 
was free to move in both radial and axial directions. Each spline centerline control point contributed two variables 
(one radial and one axial degree of freedom) to the design space. The second spline was used to vary thickness of 
the geometry about the centerline and was permitted radial motion only. Each thickness spline control point 
contributed a single variable (thickness degree of freedom) to the design space. Deformations from the centerline 
and thickness B-splines were combined to define the resultant geometry using Eqs. 12 and 13. 


R-i outer 
Ri inner 


~ ^i, 0 outer "f R r ef\- S P^^ e center (Mi center > C center') "f" s P^^ e thick(Pi thick> ^thick)\ 
~ ^i , 0 inner R r ef\- S P^^ e center (Mi center » C center ) — s P^^ e thick(Mi thick > Q/iic/c)] 


( 12 ) 


^i outer ^i inner X re f\spHne center (Mi center > ^ center )] 


(13) 


Applying Eq. 7 to Eqs. 12 and 13 yields the derivatives of the radial deformations with respect to each control point, 
denoted C center j or j. 


dRi inner 
dC center j 


dRi outer 
dC center j 


n d Spline center 

Kref ~dr — 

center j 


ui center 


( 14 ) 
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outer 
dC thick j 


— R re f 


dspline thick 

dC thick j 


ui thick 


dRi inner 
dC thick j 


— R re f 


d splivLCf-fricfc 
dC thick j 


ui thick 


(15) 

(16) 


It should be noted the derivatives in the axial direction with respect to all thickness control points ( C thick j ) are zero. 
The derivatives in the axial direction with respect to centerline control points (C center j) are given by: 

^^i inner _ ^^i outer _ y dspliue center 

^ ref -3 r> ( 1 ') 

OL centerj ° ^ center j ° L center j u i center 

At this point, the geometric parameterization has been defined using 2-dimensional polar coordinates. Extension 
to a 3 -dimensional cylindrical coordinate system representative of the actual nozzle geometry was accomplished by 
keeping the third dimensional angular coordinate, denoted 6 , constant during all deformations. Additionally, since 
the FUN3D deformation mechanics require discrete geometry specification in Cartesian coordinates, the 
transformations in Eq. 18 were applied. 


Y t = RiSin^O) 
Z t = Ricosid ) 


(18) 


Similarly, the derivatives were also transformed into Cartesian coordinates. Axial derivatives were unaffected by the 
change in coordinate system while — and — where computed using: 


dY t dY t dR t dR t 

IF = FFFF = sin ^ d t>FF 
dCj dR t dCj oCj 

dZ t dZ t dR t dR t 

dCj dR t dCj y lJ dCj 


(19) 

( 20 ) 


In Eqs. 19 and 20, is the particular partial derivative given in Eqs. 14, 15 and 16 relevant to the specific Cj 
considered. 

The FFD approach using B-splines is demonstrated for the centerbody and cowl geometries in Fig. 3. Each of 
the deformations shown is completely arbitrary to highlight method flexibility. During the actual optimization, 
deformations were only permitted downstream of the core and bypass throat areas to satisfy engine mass flow rate 
requirements during design. Constraining the entrance flow paths for the core and bypass streams also ensured 
proper geometric integration with the upstream portions of the engine. 

The complete geometry envelope for optimization is shown in Fig. 4. All surfaces were further constrained in 
parameter space to ensure zero or negative thickness regions would not occur. This was accomplished through 
careful control point placement and selection of conservative variable bounds. For example, localized thinning of the 
cowl geometry could not result in an intersection of the outer and inner cowl surfaces. Radial deformations of the 
control point influencing the centerbody tip were also not permited. In general, the cowl and shroud geometries were 
restricted from making substantial radial deformations. 
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1.5 


Center body Control 13 -Spline 



Initial and Deformed Ccnterbody Shape 





Fig 3. Arbitrary deformations to the centerbody and cowl shapes using B-spline interpolants. Perturbations to 
uniform B -splines are scaled and mapped back to deform the initial profiles. 
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Fig. 4. Baseline nozzle profile (black) can be perturbed via axisymmetic 
FFD, confined by the physical design envelope (gray). 

Example surface sensitivities to the movement of a single cowl centerline control point in the axial direction with 
respect to the x-axis are plotted in Fig. 5. Discrete analytic sensitivities were computed with respect to Cartesian x-, 
y-, and z- axes. As expected, the sensitivity appears highest in the vicinity of the control point that was perturbed. 
All geometric analytic derivatives were verified using finite difference and complex step. 



Fig. 5. Cowl surface sensitivities in the x-direction with respect to a 
single cowl centerline control point axial perturbation. 

IV. Flow and Adjoint Solvers 

The Fully Unstructured Navier- Stokes Three-Dimensional (FUN3D) 11 flow solver is the node-based CFD code 
applied in this study. FUN3D uses the finite volume method on unstructured grids to solve the discrete Euler and 
Reynolds-average Navier-Stokes (RANS) equations for compressible and incompressible flows. The solver employs 
a second-order accurate implicit point-iterative method with CFL scheduling available for convergence acceleration. 
Fluxes at cell interfaces were computed using the upwind inviscid scheme of Roe 17 . The Spalart-Allmaras 18 model 
was coupled to the mean flow equations for turbulence closure. Domain decomposition was applied to exploit high- 
performance parallel computing. 

Following convergence of the steady-state RANS equations, the set of discrete adjoint equations 19-21 were solved 
using a dual-consistent time-marching method. 22 ’ 23 One advantage of an adjoint approach is that sensitivity 
information can be obtained for an objective function with respect to many design variables at a computational cost 
comparable to that of a single flow solution. Rigorous sensitivity information from an adjoint solution may also be 
exploited to guide mesh adaptation and arrive at an, often non-intuitive, grid-converged result. This enables output- 
based engineering metrics to be captured to a high degree of precision by minimizing associated spatial 
discretization errors. 

During design optimization, volume mesh deformations were accomplished by solution of an analogous linear 
elasticity problem. 19,24 Elements near solid wall boundaries were stiffened to conform to prescribed surface 
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deflections where Young’s modulus was inversely proportional to the distance from the nearest solid boundary and 
Poisson’s ratio was assigned a uniform value of zero. All interior volume grid displacements were computed from 
the finite volume formulation of the elasticity equations using the Generalized Minimum Residual (GMRES) 25 
algorithm. 

All flow and adjoint computations were performed using 600 Intel Has well processors operating at 2.5 GHz on 
the NASA Pleiades supercomputer. Each flow or adjoint solution required approximately 1 hour of wall clock time 
to reach convergence. 


V. Optimization Algorithm 

Given the high computational cost of RANS CFD, a local gradient-based optimization approach was selected for 
efficiency. While gradient-free methods like genetic algorithms have higher probability of identifying a true global 
optimum in the presence of local minima, they require significantly more function evaluations and exhibit slow 
convergence rates detrimental for aerodynamic shape optimization. 26 

For the unconstrained optimization performed here, the standalone Sparse, Nonlinear Optimizer (SNOPT) 27 was 
implemented. SNOPT is a gradient optimizer that uses the Sequential Quadratic Programming (SQP) method. A 
Quadratic Programming (QP) subproblem is solved during each major iteration to determine the next search 
direction. The SQP algorithm implements a smooth augmented Lagrangian merit function with explicit provision for 
infeasibility in the original problem and QP subproblems. 27 The Hessian of the Lagrangian is estimated using the 
limited memory quasi-Newton method. SNOPT has been demonstrated proficient for solving large-scale 
optimization problems containing up to thousands of design variables and constraints. 

VI. Adjoint-Based Mesh Adaptation 

The baseline tetrahedral 3-D volume grid consisting of anisotropic extruded cells near all nozzle surface 
boundaries was created using the software Pointwise®. All boundary cells were resolved to a y+ < 1 with the initial 
volume containing approximately 3.5 million nodes. The z-symmetry plane grid topology and Mach number 
solution contour are shown in Figs. 6 and 7. The initial baseline grid was constructed to capture near- wall viscous 
effects only. The momentum wake, shear flow, compression and expansion features were intentionally left under- 
resolved to simplify the manual gridding process. 



Fig 6. Top: Initial baseline body-fitted grid with anisotropic extruded cells near surface 
boundaries. Bottom: Output- adapted grid to reduce the discretization error of an off-body 
pressure integral taken one cowl diameter from the nozzle centerline. 
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Baseline Grid 



Adjoint- Adapted Grid 


Fig. 7. Top: Mach number solution contour for non-adapted baseline grid. Bottom: Mach number solution contour 
for output- adapted grid. 

Starting from the baseline discretization, output-based grid refinement 28 was performed to reduce the estimated 
discretization error of an off-body pressure integral taken one cowl diameter from the nozzle centerline. This 
focused the adaptation on refining the pressure deviation from freestream in the nearfield region. High aspect ratio 
elements near all body-fitted surfaces were frozen during the refinement process to avoid disruption of the viscous 
boundary layer discretization. Additional constraints were used to control maximum anisotropy and cell count 
during refinement steps. A total of 8 refinement iterations were performed to resolve details within the plume and 
increase cell density and quality in regions impacting the nearfield pressure signature. Grid refinement operations 
included node insertion, node movement, element collapse and element swap. The final adapted grid contained 
roughly 11.4 million nodes, approximately 3.25 times the number of nodes present in the baseline grid. The z- 
symmetry plane of the adapted grid along with a corresponding Mach number flow solution is shown in Figs. 7 and 
8. While this grid was not adapted specifically to resolve thrust, the refinement increased cell densities in regions 
critical for thrust computation including the throat location of the core stream. The entire adjoint-based grid 
adaptation completed in approximately 24 hours of wall clock time utilizing 600 Intel Has well processors. 

VII. Adjoint-Based Design Optimization Results 

Starting from the adapted baseline grid, the adjoint design capabilities of FUN3D were applied to reduce the 
nearfield pressure signature one cowl diameter from the nozzle centerline. The baseline and optimized geometries 
and Mach number solution contours are shown in Figs. 8 and 9. The scale of the axes has been adjusted in Fig. 8 to 
highlight geometric differences. 

Comparison of the initial and final geometry profiles indicates the nozzle centerbody was lengthened and the 
curvature of the expansion side of the plug was reduced. The cowl was shortened and a subtle concavity was 
introduced on the outer cowl surface in contact with the freestream flow. The cowl centerline profile was adjusted 
near the exit, resulting in a higher trailing edge boattail angle for the optimized geometry. The length of the shroud 
was maintained and an increase in shroud trailing edge angle was also evident from optimization. Comparison of the 
baseline and optimized nozzle Mach number contours shows a subtle attenuation of compression features within the 
propulsive streamtube. The centerbody downstream extension resulted in diffused primary and reflected shocks 
within the core stream. The discrete flow expansion evident in the baseline solution near the cowl trailing edge is 
significantly diffused in the optimized solution. The bypass stream also appears to exit the bypass flowpath at a 
higher velocity in the optimized solution. 
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Baseline Geometry 

Optimized Geometry 



Fig 8. Comparison of baseline and optimized nozzle geometries for reducing nearfield pressure disturbances. 


Baseline Nozzle Solution 



Optimized Nozzle Solution 


Mach No. 



Fig. 9. Comparison of the baseline and optimized nozzle Mach number flow solutions for reducing nearfield 
pressure disturbances. 
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The baseline and optimized nearfield static pressure signatures taken with respect to the freestream reference 
pressure are plotted in Fig. 10. Note the optimized geometry exhibits a reduced under-pressure expansion near x = 
6.5-m. Variations to the external cowl introduced a pair of weak expansion/compression features evident between x 
= 2.5-m and 5.0-m. Given the bounds of the parameterization, it is likely these features were necessary to achieve 
the overall under-pressure reduction near x = 6.5-m. The overpressure disturbance due to the nacelle trailing edge lip 
shock which occurs near x = 7-m is also reduced when compared to the baseline signature. It is apparent the 
optimizer targeted the maximum and minimum peaks in the signature, which have the largest cost function impact. 


Baseline vs. Optimized Nearfield Pressure Signatures 


PLh 

T3 



Axial Position (m) 


Fig. 10. Comparison of the baseline and optimized nozzle nearfield pressure signatures, 
extracted one cowl diameter from the nozzle centerline. 


The thrust breakdown by component for the baseline and optimized geometries is listed in Table 2. Each value 
includes contributions of viscous forces, pressure forces and fluxes pertaining to each nozzle surface or inflow 
boundary, if applicable. As done for the geometry parameterization, the cowl and shroud were separated into inner 
and outer surfaces at the trailing edge. A 0.2% thrust increase is apparent for the optimized configuration. The 
majority of this thrust gain was due to variation in pressure forces acting on the nozzle centerbody. A small viscous 
drag increase for the plug due to axial extension was also observed, but not significant enough to offset the pressure 
force variation. 


TABLE 2: Summary of Individual Thrust Contributions (N) 

Component 

Baseline 

Geometry 

Optimized 

Geometry 

Bypass Inflow 

26478 

26473 

Core Inflow 

96228 

96238 

Inner Cowl Surface 

-1347 

-1357 

Outer Cowl Surface 

-1525 

-1552 

Inner Shroud Surface 

1505 

1536 

Outer Shroud Surface 

-1158 

-1190 

Centerbody 

-32121 

-31914 

Total 

88059 

88234 
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The optimization history is displayed in Fig. 11 and was terminated after 18 flow/adjoint cycles or 36 hours of 
wall clock time utilizing 600 Intel Haswell cores. The final nozzle design exhibited approximately a 12% reduction 
in the objective function relative to the baseline configuration. The optimizer converged predominately within just 
the first two iterations and approximately half of the control point bounds were found active upon termination. 


Nozzle Optimization History 



Fig 11. Nozzle optimization history. 


VIII. Conclusion 

Propulsion- airframe integration for ultra-low boom supersonic aircraft calls for improved methods to reduce off- 
body pressure disturbances introduced by the propulsion system or development of new methods for simultaneously 
tailoring propulsion system components in the presence of an airframe. This research aims at developing a method 
versatile enough to address both objectives. 

This study marks a first-of-its-kind effort to apply gradient-based optimization coupled with the FUN3D adjoint 
design framework to mitigate nearfield compression/expansion features for an isolated dual- stream nozzle without 
compromising nozzle performance. Careful independent tailoring of the cowl, shroud and centerbody nozzle 
surfaces was found to enable subtle reduction to compression features in the propulsive streamtube while 
simultaneously permitting a 0.2% improvement in thrust at the design point condition. The objective function, which 
was posed as a summation of near- field pressure disturbances from freestream, was reduced by approximately 12%. 

The optimization results imply that increasing the size of the design variable bounds would permit more 
geometric flexibility and further open the design space. Several challenges remain and must be addressed prior to 
performing more aggressive optimization studies. For example, permitting more substantial geometry variations 
were found to prove particularly challenging for the volume grid deformation mechanics, where negative volume 
elements were more likely to be produced. In addition, large shape changes alter the overall nozzle flow 
characteristics and deform the volume grid in a manner that numerically diffuses shock and expansion feature 
critical for nearfield signature prediction. Adjoint-based grid adaptation following each design iteration step would 
be the most desirable solution to overcome this limitation, but proves too computationally expensive at this time. 

IX. Future Work 

This research remains a work in progress and several follow-on activities are planned. Adjoint thrust derivatives 
will be implemented and validated within FUN3D. These will be used to replace the geometric bounds currently 
applied to avoid compromising core and bypass stream mass flow rates and enable geometric deformations to be 
introduced into the upstream core and bypass stream geometry. Improvements in robustness of the grid deformation 
process by solution of the linear elasticity equations will also be studied. The production of negative volume cells 
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during grid deformation frequently interrupted the design optimization process when large deformations were 
attempted. The high anisotropy of nearfield cells combined with sharp cowl/shroud trailing edges and symmetry 
plane constraints on allowable grid movement are probable causes. Alternate deformation functions may be 
considered or deformation from the baseline grid upon failure may be implemented as a work-around. The addition 
of an aggregate objective function that takes into account nozzle weight and uses volume as a potential surrogate is 
also of interest. The deformation methods developed here have been generalized and can be extended to reshape 
similar configurations, i.e. axisymmetric supersonic inlets. It is anticipated that extension to a B-spline surface-based 
parameterization will enable optimization of non-axisymmetric engine components and may also be investigated in 
future work. 
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